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Wo present the results of searches for decays of B mesons to final states with a b\ meson and a 
neutral pion or kaon. The data, collected with the BABAR detector at the Stanford Linear Accelerator 
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Center, represent 465 million BB pairs produced in e'^e annihilation. The results for the branching 
fractions are, in units of 10"^ B{B+ b+K°) = 9.6±1.7±0.9, B{B° b°iK°) = 5.1±1.8±0.5 (< 
7.8), B{B+ &+7r°) = 1.8 ± 0.9 ± 0.2 (< 3.3), and B{B° ^ bjTr") = 0.4 ± 0.8 ± 0.2 (< 1.9), with 
the assumption that B{bi ujtt) — 1. We also measure the charge asymmetry Ach{B^ bfX'^) = 
—0.03 ±0.15 ±0.02. The first error quoted is statistical, the second systematic, and the upper limits 
in parentheses indicate the 90% confidence level. 

PACS numbers: 13.25.Hw, 12.15.Hh, ll.30.Er 



Recent searches for decays of B mesons to final states 
with an axial-vector meson and a pion or kaon have re- 
vealed modes with branching fractions that are rather 
large among charmless decays: (15 — 35) x 10~^ for 
B -> ai{TT,K) [H, Q, and (7 - 11) x IQ-^ for charged 
pion and kaon in combination with a 6^ or a meson 
[1, 0] ■ In this paper we present the results of investiga- 
tions of the remaining charge states with bi accompanied 
by a tt'^ or K'^ . No previous searches for these modes have 
been reported. 

The mass and width of the bi meson are 1229.5 ± 3.2 
MeV and 142 ± 9 MeV, respectively, and the dominant 
decay is to ujir [^. In the quark model the bi is the 
/G = i± member of the J^^ = 1+", ^Pi nonet. The 
Cabibbo-favored amplitudes that mediate these decays 
are those represented by color-suppressed tree diagrams 
for the modes with tt", and "penguin" loop diagrams 
for those with K'^ . Because the 6i meson has even G- 
parity, only amplitudes in which the bi contains the 
spectator quark from the B meson are allowed, apart 
from isospin-breaking effects Q. Direct CP violation 
would be indicated by a non-zero value of the asymmetry 
Ach = (r- -r+)/(r- +r+) in the rates r±(B± ^ F±) 
for charged S-meson decays to final states F^. 

The available theoretical estimates of the branching 
fractions of B mesons to biir and biK come from calcu- 
lations based on naive factorization 0,11], and on QCD 
factorization ^ . The latter incorporate light-cone distri- 
bution amplitudes evaluated from QCD sum rules, and 
predict branching fractions in quite good agreement with 
the measurements for B 6i7r+ and B biK^ The 
expected branching fractions from QCD factorization are 
about 10x10"^ for B+ bt K° , and 3x10"^ or less for 
blK" and B ^biTT° H. 

The data for these measurements were collected with 
the BaBAR detector [13] at the PEP-II asymmetric e+e~ 
collider located at the Stanford Linear Accelerator Cen- 
ter. An integrated luminosity of 424 fb"""^, corresponding 
to (465 ± 5) X 10^ BB pairs, was produced by e+e" an- 
nihilation at the T(4S') resonance (center-of-mass energy 
^/s = 10.58 GeV). Charged particles from the e+e~ in- 
teractions are detected, and their momenta measured, 
by a combination of five layers of double-sided silicon 
microstrip detectors and a 40-layer drift chamber, both 
operating in the 1.5 T magnetic field of a superconduct- 
ing solenoid. Photons and electrons are identified with 
a CsI(Tl) electromagnetic calorimeter (EMC). Further 



charged particle identification (PID) is provided by the 
average energy loss (dE/dx) in the tracking devices and 
by an internally reflecting ring imaging Cherenkov detec- 
tor (DIRC) covering the central region. A detailed Monte 
Carlo program (MC) is used to simulate the B produc- 
tion and decay sequences, and the detector response [11|. 

The bi candidates are reconstructed through the de- 
cay sequence bi ujir, uj —> tt+tt^tt'', and 7r° 77. The 
other primary daughter of the B meson is reconstructed 
as either K'^ tt+tt" or tt*' 77. For K'^, the invariant 
mass of the pion pair is required to lie between 486 and 
510 MeV, i.e., within about 3.5 standard deviations of 
the nominal K'^ mass [s!] . The minimum energy for a tt"- 
daughter photon is 30 MeV (50 MeV for a primary tt"), 
and the minimum energy of a tt*^ is 250 MeV. The in- 
variant mass of the photon pair is required to lie between 
120 and 150 MeV, or within about two standard devia- 
tions of the nominal tt*^ mass. For the bi and to, whose 
masses are treated as observables in the maximum like- 
lihood (ML) fit described below, we accept a range that 
includes wider sidebands (see Fig.[T|). Secondary charged 
pions in 61 and to candidates are rejected if classified as 
protons, kaons, or electrons by their DIRC, dE/dx, and 
EMC PID signatures. For a candidate we require a 
successful fit of the decay vertex with the flight direc- 
tion constrained to the pion pair momentum direction, 
that yields a flight length greater than three times its 
uncertainty. 

We reconstruct the i?-meson candidate by combin- 
ing the four-momenta of a pair of primary daughter 
mesons, using a fit that constrains all particles to a 
common vertex and the tt*^ mass to its nominal value. 
From the kinematics of T{AS) decay we determine the 

energy-substituted mass niEs — \J \s — and energy 

difference IS.E = Eb — \^/s^ where {Eb,Pb) is the B- 
meson four-momentum vector, and all values are ex- 
pressed in the T(45') rest frame. The resolution in tties is 
2.4-2.8 MeV and in AE is 22-46 MeV, depending on the 
decay mode. We require 5.25 GeV < mES < 5.29 GeV 
and \AE\ < 100 MeV. 

We also impose restrictions on the helicity-frame decay 
angles of the bi and uj mesons. The helicity frame of a 
meson is defined as the rest frame of the meson with 
z axis along the direction of boost to that frame from 
the parent rest frame. For the decay bi ujtt, Oh^ is 
the polar angle of the daughter pion, and for lo — > 37r, 
6,^ is polar angle of the normal to the 37r decay plane. 
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Since many misreconstructed candidates accumulate in a 
corner of the cos^b^ vs cosOi^ plane, we require cos^fc^ < 
min(1.0, 1.1-0.5 x Icose*^]). 

Backgrounds arise primarily from random combina- 
tions of particles in continuum e+e" qq events {q — 
u, d, s, c). We reduce these with a requirement on the an- 
gle 9t between the thrust axis I4I of the B candidate in 
the T{AS) frame and that of the charged tracks and neu- 
tral calorimeter clusters in the rest of the event (ROE). 
The event is required to contain at least one charged track 
not associated with the B candidate. The distribution is 
sharply peaked near |cos6't| — 1 for qq jet pairs, and 
nearly uniform for i3-meson decays. The requirement, 
which optimizes the expected signal yield relative to its 
background-dominated statistical error, is | cos^tI < 0.7. 

The average number of candidates found per event in 
the selected sample is in the range 1.3 to 1.6 (1.4 to 1.6 
in signal MC), depending on the final state. We choose 
the candidate with the largest confidence level for the 
S-meson geometric fit. 

In the ML fit we discriminate further against qq back- 
ground with a Fisher discriminant that combines five 
variables: the polar angles, with respect to the beam 
axis in the T(4S') rest frame, of the B candidate momen- 
tum and of the B thrust axis; the flavor tagging cate- 
gory; and the zeroth and second angular moments Lq^2 
of the energy flow, excluding the B candidate, about the 
B thrust axis. The tagging category [l^] is the class 
of candidate partially reconstructed from the ROE, de- 
signed to determine whether, in a signal event, it rep- 
resents a B OY B meson. The moments are defined by 
Lj — ^j^Pi X jcos^il"' , where 9i is the angle with respect 
to the B thrust axis of track or neutral cluster i, Pi is 
its momentum, and the sum excludes the B candidate 
daughters. The Fisher variable provides about one stan- 
dard deviation of separation between B decay events and 
combinatorial background. 

We obtain yields for each channel from an extended 
ML fit with the input observables Ai5, toes, -^j and the 
resonance masses rub^ and . The selected data sample 
sizes are given in Table HI Besides the signal events these 
samples contain qq (dominant) and BB with b ^ c 
combinatorial background, and a fraction of cross feed 
from other charmless BB modes, which we estimate from 
the simulation to be (0.5-1.1)%. The last include non- 
resonant WTTTT, LoK-K, and modes that have final states 
different from the signal, but with similar kinematics so 
that broad peaks near those of the signal appear in some 
observables. We account for these with a separate com- 
ponent in the probability density function (PDF). 

The likelihood function is 



N 



c = cxp -^y, HE^^ 



(1) 



V, (mEslV, {r)V, iAE^)V, {ml )V, (to^). 



where N is the number of events in the sample, and for 
each component j (signal, combinatorial background, or 
charmless BB cross feed) , Yj is the yield of events and 
Vj{x^) the PDF for observable x in event i. The sig- 
nal component is further separated into two components 
(with proportions fixed in the fit for each mode) repre- 
senting the correctly and incorrectly reconstructed can- 
didates in events with true signal, as determined with 
MC. The fraction of misreconstructed candidates is 32- 
40%, depending on the mode. The factored form of the 
PDF indicated in Eq. [l]is a good approximation, partic- 
ularly for the combinatorial qq component, since we find 
correlations among observables in the data (which are 
mostly qq background) are generally less than 2%, with 
none exceeding 5%. The effects of this approximation 
are determined in simulation and included in the bias 
corrections and systematic errors discussed below. 

We determine the PDFs for the signal and BB back- 
ground components from fits to MC samples. We cal- 
ibrate the resolutions in AE' and toes with large data 
control samples of B decays to charmed final states of 
similar topology (e.g. B — > D{KTnT)TT, B — > D{Knn)p). 
We develop PDFs for the combinatorial background with 
fits to the data from which the signal region (5.27 GeV < 
TOES < 5.29 GeV and |A£;| < 75 MeV) has been ex- 
cluded. 

The functions Vj are constructed as linear combina- 
tions of Gaussian and polynomial functions, or in the 
case of Toes for qq background, the threshold func- 
tion xvT^-^exp [— ^(1 — x^)] , with argument x = 
'^m-Es/y/s and shape parameter ^. These functions are 
discussed in more detail in [l5|, and are illustrated in 
Figs. [1] and H 

We allow the parameters most important for the deter- 
mination of the combinatorial background PDFs to vary 
in the fit, along with the yields for all components, and 
the signal and qq background asymmetries. Specifically, 
the free background parameters are: ^ for toeSj linear 
and quadratic coefficients for AE, and the mean, width, 
width difference, and polynomial fraction parameters for 

We validate the fitting procedure by applying it to 
ensembles of simulated experiments with the qq compo- 
nent drawn from the PDF, into which we have embedded 
known numbers of signal and BB background events ran- 
domly extracted from the fully simulated MC samples. 
By tuning the number of embedded events until the fit 
reproduces the yields found in the data, we determine 
the biases that are reported, along with the signal yields, 
in Table U 

In Figs. [1] and [2] we show the projections of the PDF 
and data for each fit. The data plotted are subsamplcs 
enriched in signal with the requirement of a minimum 
value of the ratio of signal to total likelihood (computed 
without the plotted variable) that retains (30-50)% of 
the signal, depending on the mode. 
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TABLE I: Number of events A'' in the sample, fitted signal yield Ys, and measured bias (to be subtracted from Ys) in events 
(ev.), detection efficiency times secondary decay branching fractions e, significance S (with systematic uncertainties included), 
and branching fraction and charge asymmetry with statistical and systematic error. 



Mode N (ev.) Ys (ev.) Bias (ev.) e (%) S (a) B (IQ-^) A, 



ch 



btK° 9841 164tll 15±7 3.4 6.3 9.6 ±1.7 ±0.9 -0.03 ± 0.15 ± 0.02 

6?A'° 5420 58t{r 5±3 2.2 3.4 5.1 ± 1.8 ± 0.5 (< 7.8) 

btn° 28787 7lig 8±4 7.7 1.6 1.8 ± 0.9 ± 0.2 (< 3.3) 

fe?7r° 10554 6ll6 -2±2 4.8 0.5 0.4 ± 0.8 ± 0.2 (< 1.9) 




FIG. 1: Distributions for signal-enhanced subsets (see text) 
of the data projected onto the fit observables for the decay 
B+ b+K"; (a) thes, (b) AE, (c) JF, (d) m(7r+7r-7r'') for 
the uj candidate, and (e) m{(jjn) for the bi candidate. The 
solid lines represent the results of the fits, and the dashed 
and dot-dashed lines the signal and background contributions 
respectively. 



FIG. 2: Distributions for signal-enhanced subsets (see text) of 
the data projected onto niES (a, c, e) and AE (b, d, f) for the 
decays B° b^K" (a, b), B+ b+Ti" (c, d), and B° ^ 6^° 
(e, f). The solid lines represent the results of the fits, and 
the dashed and dot-dashed lines the signal and background 
contributions respectively. 



We compute the branching fraction by subtracting the 
fit bias from the measured yield, and dividing the result 
by the number of produced BB pairs and by the efficiency 
times B{uj vr+Tr-Tr") = 89.1 ±0.7% (and for the modes 
with K'^, B{K° K^s ^ TT+TT") = i(69.20 ± 0.05)%) 

. The efficiency is obtained from the MC signal model. 
We assume that the branching fractions of the T{AS) to 
B^B^ and B^B^ are each equal to 0.5, consistent with 



measurements Q- The results are given in Table U along 
with the significance, computed as the square root of the 
difference between the value of —2 In £ (with additive sys- 
tematic uncertainties included) for zero signal and the 
value at its minimum. 

Systematic uncertainties on the branching fractions 
arise from the PDFs, BB backgrounds, fit bias, and effi- 
ciency. PDF uncertainties not already accounted for by 
free parameters in the fit are estimated from the consis- 
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tency of fits to MC and data in control modes. Varying 
the signal-PDF parameters within these errors, we esti- 
mate yield uncertainties of (1.6-6.4)%, depending on the 
mode. We estimate the uncertainty of the MC model 
of misreconstructed signal by performing alternate fits 
with a signal PDF determined from true signal events 
only; we find differences of 1-4 events between these and 
the nominal fits. The uncertainty from fit bias (Table 
HI includes its statistical uncertainty from the simulated 
experiments, and half of the correction itself, added in 
quadrature. For the BB backgrounds we vary the fixed 
fit component by 100% and include in quadrature a term 
derived from MC studies of the inclusion of a 6 ^ c com- 
ponent with the dominant qq background. Uncertainties 
in our knowledge of the efficiency include 0.5% x Nt and 
1.5% X N^, where Nt and are the numbers of tracks 
and photons, respectively, in the B candidate. The un- 
certainties in the efficiency from the event selection are 
below 0.5%. 

We study asymmetries from the track reconstruction 
(found to be negligible), and from imperfect modeling of 
the interactions with material in the detector, by mea- 
suring the asymmetries in the qq background in the data 
and control samples mentioned previously, in comparison 
with MC [l3| . We assign a systematic error for Ach equal 
to 0.01. 

With the assumption that B{bi ujtt) — 1, we obtain 
for the branching fractions (in units of 10^^): 



B{B+ - 




= 9.6 ±1.7 ±0.9 


B{B"- 


^ fe?i^") 


= 5.1 ±1.8 ±0.5 (< 7.8) 


B{B+ - 




= 1.8 ±0.9 ±0.2 (< 3.3) 


S(B°- 




= 0.4 ±0.8 ±0.2 (< 1.9) 



The first error quoted is statistical and the second sys- 
tematic. We find no evidence for the modes with tt*^; 
the evidence for B{B^ biK'^) has a significance of 
3.4 standard deviations. For these modes we quote also 
90% confidence level upper limits, given in parentheses. 
We observe the decay B{B^ biK*^), and measure the 
charge asymmetry 

Ach{B+ ^bfK") = -0.03 ±0.15 ±0.02. 

The QCD factorization estimates for the branch- 
ing fractions and charge asymmetry (0.014) agree with 
these measurements within experimental and theoretical 
errors. We find no evidence for direct CP violation in 
B{B+ -> 6+i^"). 
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